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Abstract 



I Decoherence in quantum systems which are classically chaotic is studied. It is well- 



known that a classically chaotic system when quantized loses many prominent chaotic 
traits. We show that interaction of the quantum system with an environment can 
I under general circumstances quickly diminish quantum coherence and reenact some 

characteristic classical chaotic behavior. We use the Feynman- Vernon influence func- 
tional formalism to study the effect of an ohmic environment at high temperature on 
two classically-chaotic systems: The linear Arnold cat map (QCM) and the nonlinear 
Q^' quantum kicked rotor (QKR). Features of quantum chaos such as recurrence in QCM 

^ • and diffusion suppression leading to localization in QKR are destroyed in a short time 

due to environment-induced decoherence. Decoherence also undermines localization 
and induces an apparent transition from reversible to irreversible dynamics in quan- 
^ ■ tum chaotic systems. 



1 Decoherence and Disappearance of Recurrence in 
the Quantum Cat Map 

Arnold's cat map is a linear area-preserving map T on a torus in phase space formed by 
identifying the boundaries of the interval [0, 2n] in both the coordinate Q and the momentum 
P directions |^. (Because of this the area of the torus is characterized by Planck's constant 
which takes on the values h = 27r/A/', where J\f is the number of sites in both the coordinate 
and the momentum directions in the phase space.) From time step j to j ' + 1 it is given by 

Q^+i\_(^ b \ ( QA -t( (1.1) 



Pj+i / \ c d I \ Pj I \ Pj 
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where detT = 1 guarantees area preservation. The degree of chaos depends on the choice 
of T. The eigenvalues of T are either both real or both imaginary. In the latter case, T is 
elliptic, the motion becomes periodic and no sensitive dependence on the initial condition is 
observed. When T is hyperbolic, the motion is chaotic. 

Quantized cat map is studied in detail by Hannay and Berry 0. The matrix has to 
assume a special form in order to yield nontrivial values of the progagator for the map. We 
choose 

for the elliptic and the hyperbolic cases respectively. 

For the special choice of the matrix elements Ti, T2 made above the propagator takes on 
the simple forms, 

Ui{j + l,j) = ^exp[-^Q,Q,.+i], f/2(j + l,j) = ^exp[^(Q|-Q,Q,+i + Q2^J]. (1.3) 

Since each iteration describes a permutation among sites, each site belongs to a periodic 
orbit. Thus the quantum dynamics follows the classical way, resulting in the recurrence of 
the wave function (or equivalently, the Wigner function |^). 

We now couple the system linearly to a bath of harmonic oscillators with coordinates 
qa and momentum {a = 1, ..N) described by the Hamiltonian Hb and the interaction 
Hamiltonian He 

^b = E(t + ^)' ^c = Y.C.Qcio.. (1.4) 

a=l ^ «=1 

where Q is the coordinate of the system and Ca is the coupling constant of Q to the ath 
oscillator in the bath. By integrating out the bath variables, we get the reduced density 
matrix, 

Pr(g„g;-,t) = / ntWMg;exp^[5(g)+5c(g,g.)+5ij(ga)-5(g')-^c(g',g:.)-5B(g;)]. 

(1.5) 

where S is the classical action of the system which appears as the exponent of the propagator 
in (1.3). Sb-, and Sq are the actions for bath and interaction, respectively. The propagator 
Jr for the reduced density matrix from time steps j to j + 1 is 

^r(g,+i,g;+i i g„g;,t) = / dqdq' ^^^'-\s{q) - s{q')^a{q,q% (i.e) 

in a path-integral representation |P, ^, ||, where 

^A(g,g') = -^ t ds rds'r{s)[-zfiis-s')R{s')-u{s-s')r{s')] (1.7) 
h h Jo Jo 



is the influence action. Here r = ,R = and /i(s),z/(s) are the dissipation and 

noise kernels respectively 0]. 
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If we consider the simplest case of an ohmic bath at high temperature kT > fiA » fiooa 
0], and consider times shorter than the relaxation time, then we obtain a Gaussian form 
for the influence functional, with ^A{Q, Q') = — ^^^1^'^ ^j^f- where the noise kernel becomes 
local = 2M'-fkT6{s) and 7 is the damping coefficient. The unit-time propagator becomes 

j,(g,+i,g;.+i | g,-,g;) = (j,(g,+i,g;+i | g,-,g;,o) = (exp^[5(g,+i,g,)-5(g;+„g;)+eri+i]). 

(1.8) 

Here ^ is a Gaussian white noise given by 



(0 = 0, (exp^er)=exp[-^^^r2] (1.9) 



where ( ) denotes statistical average over noise realization ^. 
For the elliptic map, we get 

J.(g,+i,g;.+i I g,-,g;.,0 = (-^)^/'exp[^(-r,i?,+i -r,+ii?, +er,+i)]. (1.10) 
and for the hyperbolic map, 

MQj+i, Q'j+i I Qj, Q'j,0 = i^y^'^ exp[^(2rji?j + 2rj+iRj+i - rjEj+i - rj+iRj + ^rj+i)] 

(1.11) 

The Wigner function is defined as 

1 2i 
W(R,p) = — ^(R + rU*(R-r)exp(—pr)dr. (1.12) 
nn Jo n 

where p is the momentum conjugate to r. The propagator K for the Wigner function is 

KiRj+i,Pj+i I Rj,Pj,0 = K,K,+Jr{Qj+i,Q'j+i I Qj,Q'j,0(^wf{Pjrj-Pj+irj+i). 

(1.13) 

This is reduced to the form of the classical cat map. For the elliptic case, 

= -Pj+i + Pj = Rj+i- (1-14) 

For the hyperbolic case, 

Rj = -Pj+i + 2i?,+i + Pj = -3i?,+i + 2pj - 2e (1.15) 

Without noise, quantum evolution follows classical permutation the phase space is 
divided by a finite number of different periodic orbits and the period is known to increase 
roughly proportional to Af with some irregular oscillation. When coupled to a bath, the 
cat map is exposed to a Gaussian noise in each time step. The discretized noise induces 
transitions between different periodic orbits in an irregular way. Interaction with an envi- 
ronment blurs the recurrence of physical quantities in the quantum map. Fig.l shows Trp^, 
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the linearized entropy (with the reversed sign) for various cases. If there is no interaction 
with the environment, the entropy is constant for both regular and chaotic cases. Quantum 
recurrence is evident even when the system is chaotic. When interaction sets in, Trp^ decays 
exponentially, showing that the system rapidly decoheres. The rate of decoherence is much 
faster in chaotic systems than in regular systems ||^. It suggests that recurrence would be 
less evident in a decohering chaotic system. In Fig. 2, we show the mean displacement of 
points in the phase space as a function of time steps. This is defined by / = ^ (Ax^ + Ap^), 
where Ax and Ap are the displacements from the initial phase space points, and ( ) denotes 
averaging over noise distributions. In the chaotic case, we see that recurrence disappears 
with just a small amount of noise (Fig. 2a) whereas in the regular case, the same amount 
of noise does not alter the qualitative picture of recurrence (Fig. 2b). In both cases, the 
decohered quantum system behaves close to the classical picture in which the regular and 
chaotic dynamics are clearly distinguished. In spite of the discreteness of the points on the 
torus, the system behaves effectively classically due to the influence of the environment. 



2 Decoherence and Delocalization in the Quantum Kicked 
Rotor 

The kicked rotor is one of the most intensively studied models from both the quantum and 
classical point of view ||^. The Hamiltonian of the kicked rotor is given by 

2 

H=^ + ^cosxS°^_^5(t - j) (2.1) 

which describes a one- dimensional rotor subjected to a delta-functional periodic kick at 
t = j. Here x is the angle of the rotor with period 27r, m is the moment of inertia, p is 
the angular momentum, and K is the strength of the kick which measures the nonlinearity. 
When K > 1, the system becomes chaotic over the entire phase space. 

The quantum dynamics of the kicked rotor is given by the corresponding Schrodinger 
equation 

d 

ih—^{x,t) = - — —iP{x,t) + K cos xE,6{t - j)^{x,t) (2.2) 

where ip is the wave function of the rotor. 

Denoting ipj as the wave function ijj{x,t) at each discrete time t = j, and integrating 
(2.2) from j to j ' + 1, we obtain 



; / \ r ■ ^ ^ r .-^COSX, , , , ,„ 

tj+i{x} = ^M-'^^-Q^l exp[-2 — \ipj{x} (2.3) 



The quantum kicked rotor (QKR) is known to exhibit dynamical localization. After some 
relaxation time scale, the wave function becomes exponentially localized in the momentum 
space 13] . This may be interpreted as a particle moving in a lattice with a quasi-random 
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potential. This heuristic picture seems to justify the analogy between the quantum kicked 
rotor to the tight binding model with an exponentially decaying hopping parameter which 
is known to show Anderson localization [Q. Dynamical localization in this context arises 
from the suppression of classical diffusive behavior by the quantum dynamics. However, 



as shown by Ott, et.al. |Ty], a small external noise can break the localization. Sufficient 
amount of noise would induce the quantum system to exhibit classical diffusive behavior. 
Dittrich and Graham studied this problem |jTl| by coupling the system to a zero temperature 



harmonic oscillator bath and analysed solutions to the master equation. Cohen and Fishman 



presented the most detailed study of this problem for an ohmic bath |T2[. Here we want 
to approach these issues from an environment-induced decoherence point of view [|13|. We 
begin by calculating the density matrix for the kicked rotor coupled to an environment. 

We introduce a linear coupling of the system momentum p with each oscillator coordinate 
qa{(y = 1, ..A^) in the bath in the form Hq = ^a=i^a<laP (Here q,p without the subscript 
a denote the system coordinate and momentum variables). As before, we assume an ohmic 
bath and examine the time period where dissipation is small. Under these assumptions, the 
unit time propagator for the wave function U^{j + 1, j) is given by 

2 

U({j + 1, j) = 6^P[~^^] exp[--ircosx] exp[--^p] (2.4) 

where, as before, the noise term ^ arises from using a Gaussian identity in the integral 
transform of the term involving the noise kernel in the influence functional. Summing over 
all noise realizations ( ) gives the desired reduced density matrix, 

PrjiP, P') = ii^^ip) ^j,^{.P'))^ (2.5) 

where 

V',+i,^(p) = f/5(j + l, J)^,,^(P) (2.6) 

Loss of quantum coherence is measured by the density matrix becoming approximately 
diagonal. Trpl can be expressed as 

Trpl = (S,S,,^5(p) r^ip') i;^ip') r^ip))^,^' (2.7) 

where ( )^^^' denotes the statistical average of all possible noise histories of two independent 
noises ^, ^' defined at each time interval from j to j + 1. At high temperatures ^(t), ^'(t) are 
reduced to two time-uncorrelated independent Gaussian white noises defined at each time 
step. 

We see that there is a close relation between the breaking of dynamical localization 
and quantum decoherence. In Fig. 3 we plot the linearized entropy Trp^ versus the energy 
(p^). This shows that delocalization occurs as quantum coherence breaks down, suggesting 
that delocalization and decoherence occurs by the same mechanism. As the nonlinearity 
parameter K increases, the system decoheres more rapidly. At the same time, the amount 
of delocalization measured by the diffusion constant increases. 
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This may be explained in the following way: Because the coupling is through the mo- 
mentum, the noise term does not involve any nonlinearity. The time scale for the system 
to lose coherence is given hj tr, = {XtdB / ^vY / 1 ^ where XtdB = h/\/2TTmkT is the thermal 
de Broglie wavelength, and 5p is the relevant momentum scale. After this time, noise will 
destroy the quantum coherence between such momentum separations. In the kicked rotor 
case, localization will occur due to the coherence around 5p ~ A, where A ~ //l is the local- 
ization length. Since I ~ K^, this gives t_D ~ K~'^. This shows that nonlinearity increases 
the rate of decoherence. 



The relation between the diffusion constant D and the noise strength is given in ||T0|, [12 
For our case, K/h ^ 1 and for weak noises, we can consider the particle as undergoing a 
random walk with hopping parameter 1/tc- Then D = A'^/tn = (A^/Aj^^)/7. 



3 Decoherence and Irreversibility in Quantum Chaos 

The Wigner function is often used to examine the quantum to classical transition. [] The 
Wigner function at time t = j is defined as 

+2n 



WjiX,p) = ^ I dye^^y p,{x+y-,x-y-\ (3.1) 



where X = |(x + x'),y = x — x' . From (2.3), the unit-time propagator for the Wigner 
function of the QKR is found to be 

W,+,{X,p) = e-'^^'^e-P^-W,{X,p) (3.2) 

where Ap = e^^^ — e~^^^ measures the effect of the kick. We can see the effects of quantum 
corrections is seen more clearly if we expand Ap in orders of h: 

The first exponential contains the classical propagator and the second contains quantum 
corrections of even orders of h. Thus we get 

Wj+i{X,p) ^ e^^''''''^'^Wj{X -{p + Ksinx),p- Ksinx) (3.4) 

where the Wigner function with the new arguments depicts classical evolution. This map 
alone is the source of streching and folding of volume in phase space which signify classical 
chaos. 



^For a linear system the Wigner function is known to show a smooth convergence to the classical Liouville 
distribution. But if the system Hamiltonian has a nonlinear term, quantum corrections associated with the 
higher derivatives of the potential pick up the rapid oscillations in the Wigner function and it no longer 
has a smooth classical limit However, upon interaction with an environment, a coarse-grained Wigner 
function can have have a smooth classical limit p5[ for nonlinear systems. 
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If the initial system wavefunction is described by a Gaussian wave packet with width 
Sp{>> h), we would expect to see a classical-like evolution of the packet at short times. 
When the width of the contracting wave packet gets so small as comparable to h, the effect 
of quantum corrections from higher h order terms in (3.4) set in. By comparing the classical 
and quantum terms, we see that quantum corrections will become important when 6p{t) ~ h. 
Here 6p{t) = (5p(0)e~^*, where the Lyapunov exponent A ~ ln{K/2) . Thus we can deduce 
the Ehrenfest time f\ for QKR to be ~ j^i^^^^ Note that in the continuum case, this 



definition gives us a different time scale for each term in the expansion [|T6 . 

The major effect of the bath (at times short compared with the relaxation time) is the 
appearance of a diffusion term in (3.4), 

Wj+i{X,p) ^ e^^^e^^''''^^piy,(X -{p- Ksinx),p- Ksinx) (3.5) 

Competition amongst the three terms with different physical origins is apparent: The first 
term in (3.5) is the quantum diffusion term, the second is the quantum correction term, and 



the third is purely classical evolution. As discussed by Zurek and Paz |]T6|, if D is sufficiently 
large, the effect of quantum corrections becomes inconspicuous. In this case, the diffusion 
term traces out a small scale oscillating behavior before quantum corrections have a chance 
to change classical evolution. Then one may expect the time evolution of the Wigner function 
to be like that of classical evolution with noise. The role of quantum diffusion is to add some 
Gaussian averaging so that the contracting direction in phase space will be suppressed while 
it does not affect the stretching direction. As long as the width of the wave packet is large 
such that the first term is negligible, the evolution should be Liouvillian (time reversible if 
we assume infinite measurement precision). Furthermore, we expect that after the width of 
the packet along the contracting direction becomes comparable to the diffusion generated 
width (in the Gaussian wave packet) , the dynamics will start showing irreversible behavior 
arising from coarse graining (as distict from irrreversibility from instability). Consequently, 
entropy should increase in this regime. In Fig. 4a, we plot the von Neumann entropy for 
the dynamics of (3.5). We can see three qualitatively different regimes: I. the Liouville 
regime: the entropy is constant and the dynamics is time reversible. II. the decohering 
regime: the entropy keeps increasing due to coarse graining. III. the finite size regime: due 
to the bounded nature of the phase space, the entropy shows saturation. Our result from 
quantitative analysis seems to confirm the qualitative description of Zurek and Paz [^] who 
used the inverted harmonic oscillator potential as a generic source of instability. Since the 
phase space in their model is not bounded they do not see Regime III. Similar features 
appear in the quantum cap map (Fig. 4b) In this case, the full quantum dynamics can be 
calculated in a simple way. Resemblance with the result of a classical rotor with noise is 
obvious. However, in this case, the stable entropy is smaller than the maximum value which 
may be explained as a finite (phase space) size effect. 



^Ehrenfest time is customarily defined as the time when quantum dyanmics can be adequately described 
by the classical equations, i.e., the time t < when the Wigner function or the expectation value of any 
observable follow classical trajectories. 
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Figure Captions 

Figure 1 The hnearized entropy (with reversed sign) Trpl is plotted here as a function of 
time. If there is no environment, the entropy is constant for both hyperbolic and elliptic 
cases, indicating the purity of the state. For the hyperbolic map, even though classically 
this system is strongly chaotic, the corresponding quantum system does not show chaotic 
behavior. This situation changes drastically when the system interacts with a thermal bath: 
entropy keeps increasing due to coarse graining. Note that in the hyperbolic case (solid line) 
the rate of entropy increase is greater than in the eUiptic case (dotted line) . N = 50 is used 
here (also in Fig.2). 

Figure 2 The mean phase space point displacement is shown. When there is no environment 
(dotted line), the system shows recurrence in both hyperbolic (a) and elliptic (b) cases. In 
the presence of an environment, the hyperbolic map loses the recurrence behavior (solid line) 
under a Gaussian noise with a — 0.08 and maintains a near-constant value, indicating the 
ergodicity of the classical map. On the other hand, the ellptic map still shows recurrence 
with the same amount of noise, suggesting classical periodicity. 

Figure 3 Trpl (solid line, left scale) and < > (dashed line, right scale) are plotted against 
time for K — 12 and Ti — 1.52. The upper solid hne and the lower dashed hne correspond to 
the case when there is noise, with a — 0.5. As the noise strength increases to cr = 1.5, the 
decoherence time shortens, and Trpi decays rapidly (the lower solid line). This accompanies 
the increase of diffusive behavior in < > (upper dashed line). 

Figure 4 The von Neumann entropy is plotted versus time for (a) the quantum kicked rotor 
with an environment. Here, h — 1.52, a — 0.08 and K — 1.2. Entropy stays at zero (re- 
versible dynamics) until a transition regime, after which the dynamics becomes irreversible, 
(b) the quantum cat map, with the same parameters and the same amount of noise. We see 
the same qualitative feature as in the QKR case. 
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